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We present difFerential equation for evolution of interface based on continuous approximation of 
Temkin's model. 

We start from a free energy functional of type^: 

F[<j>] ^ yaV01n(aV0) + (/.2/^ + (l-<^)'/'^ (1) 
+ 20(1 - 

Where (j> is the order parameter equivalent to fraction of solid, and 1 — </) is the fraction of liquid in Temkin's model 
of solidification^, we can also assume that (f> represents fraction of phase A and 1 — (f> fraction of phase B in an AB 
system. Consequently f^, and f^^ are free energies of pure solid, liquid and liquid-solid phases, respectively. In 
case of Temkin's model a is vector representing crystal lattice parameter, but in above functional we can generalize 
this parameter to represents a characteristic length of growth not necessarily related to the lattice parameter. The 
lattice vector a can be some function of the order parameter which controls the evolution of the characteristic length 
and depend on additional space variables including Euler angles but this further complicates our model. By doing 
the variational derivative of Eq. ^ and defining a/A = f^^ — ^ ^-^ and P/2 — f^^ — we have: 

X {aldl + a^dl + aldl + 2alald,y + 2alald^, + ^alaldy^) - + = (2) 

Eq. ^ describes the general case of anisotropy growth of interface in three dimensional space. The parameters a 
and /3 are defined as L/kT and Aji/kT, respectively. The Jackson parameter a is measure of surface energy and (3 is 
difference in chemical potential between two phases^. 

We restrict our considerations only on isotropic case where: ax — ay — az — a and one dimensional case reducing 
Eq.dSl) to: 

We will now give some solutions of previous evolution equation. 

I. EQUILIBRIUM SOLUTION 

The solution of differential equation ^ for equilibrium conditions /? = and ^ = is: 4>{x) = — tanh(a/2a(a; — (5)) 
where 5 is an arbitrary constant representing current position of the interface. The thickness of the interface is given 
by: 



2a 

a 



(4) 



This result is in agreement with Temkin's theory and it's already given in^. It predicts that interface thickness would 
depend on the parameter which Temkin defines as: 



ze L 
kf ~ kf 



(5) 



Where z is number of nearest neighbors or coordination number and e — egi — is binding energy of the interface. 

Interesting consequence of relation ([5|) is prediction of the thinner interface film for higher energy of binding when 
the thickness d is normalized with lattice constant a. This can be compared with experimental facts for different 
materials. 

Using Eq.([5]) we can estimate the equilibrium thickness of the intergranular film for usual ceramic materials. If 
we take approximate values for ceramics materials, the equilibrium temperature of about 2200 K and corresponding 
value for latent heat L = 1.5 * 10"^" J/atom we get the thickness of intergranular film which is in agreement with the 
experimental findings and previous theoretical work of Bobeth et ali^. 
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II. VECTOR CASE 



Providing that we have slow growth of the crystal in x-direction we expect a slow change of Qx due the small 
temperature gradient. Thus in the limiting case of the adiabatic change we can change with the constant value of 
a. In the general case will change its value with the time ax{t). We would like to describe the change of as a 
function of time and the simplest model we can imagine is a harmonic oscillator. Then we can define as: 

— a(l + cos(ti)t)) (6) 

Where lo is the mean angular velocity of all atoms at the interface. Solution of our new differential equation in the 
case of vanishing chemical potential /3 = is: 

(j){x,t) — — teLnh{m{x — 6)) X (7) 
2am{a^ + uj"^ + a'^ cos{ujt) + aoj sin{ujt)) 

Where m ^ N Eq.(I7]) comes from the solution by separation of variables. This is wave-like tangent hyperbolical 
solution contains two parts, first which is proportional to thickness of the interface and second is harmonic oscillations 
of atoms at the interface. The harmonic oscillations of the diffuse interface will be govern by ratio a/u. 



III. EQUILIBRIUM SOLUTION IN 2 D CASE 

In this section we map the solution of differential Eq. ^ in the case of cylindrical symmetry to solution of Burgers 
vector equation: 

^ + UVU = lyV'U (8) 
ot 

This non-linear partial differential equation is described and solved in the detail Main idea is to change the 

vector field U with the gradient of the scalar field -0 imposing V x [/ = and VV-" = U. Modifying Eq. ^ we can 
simplify equation to: 

f + (^)^-VV = i^W (9) 
Where E{t) is function of time only. Applying the same trick to our Eq. ^ we get simplified version: 

aV2u-|(Vu) + /3Vu-^Vu = (10) 

Where u = \7(j) and V x m = 0. For case of the cylindrical symmetry and case where lattice vector a is always parallel 
to normal of phase field, Eq. (fTUl) can be reduced to: 

„M|w))_<.,,,_^ ,11) 

dr 

Equilibrium solution to Eq. (fTTj) is given as: 

2a/a 



(12) 



[r — S) ln(r — S) 

As in the one dimensional case the thickness is characterized by 2a/ a, but it has logarithmic divergence to infinity. 



IV. DEPENDENCE OF THE INTERGRANULAR FILM THICKNESS ON THE MISFIT ANGLE AND 

GRAIN BOUNDARY ANGLE 

The potential energy of grain boundary in Read-Shockley only depends on the dislocations in the crystal, produced 
by the mismatch of the crystal lattices. If the amorphous intergranular film is in the thermodynamic equilibrium with 
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crystalline solid, the necessary condition for thermodynamic equilibrium is the minimum of the free energy leading 
to conclusion that amorphous intergranular film has a lower free energy than grain boundary between two crystalline 
solids. If we take that is angle of grain boundary and 6 is misfit angle, then from Read-Shockley theory of the grain 
boundary dislocations^ in the crystals we can write potential energy of grain boundary as: 

E = EoB{A-\ne) (13) 

Necessary condition for existence of the amorphous intergranular layer is: 

E ^ ^amorph (-^^) 

Where: 

Eo = G'a(cos(^ + sin<y£')/(47r(l — c)) 
Eamorph —energy of an amorphous-crystal interface 
A —function depending on (p 
<f> —misfit angle 

We assume that energy of amorphous layer is exactly equal to grain boundary energy of the Read-Schokley theory, 
which fulfil the necessary condition for thermodynamic equilibrium of amorphous film. In perfect solid the lattice 
vector a is oriented in one particular direction, while in the diffuse interface the lattice vector has statistical distribution 
. This can be pictured as Ising XY-model where grain boundaries are regions in the plane where spins have different 
orientations and regions of the crystal are perfectly oriented. 

The probability of an orientation of a lattice vector in the space angle is given by: 

dWo = const dK (15) 

Or using the Boltzmann factor: 

dW = const e^"^^^^^''^^^'''^ sinipdip (16) 



/ a^dW 

-Jdw~ 

J a COS if exp[^^i^|^^^^pi0(A - In 9) /kT] sin ipdip 

I exp[^^ii(|^^i^0(A - ln0)/kT] sin^d^ 
aV2/2 



If we assume that cos ip ^ simp which is true for small angles, the exponent in exponential function is simply 
Ga{cos(p + sini^) « Gacosi^, changing the integration variable x = a cos ipG/ An {1 — (t)9{A — \n0) and collecting all 
parameters in A = 6{A — In 6*) we have the following integral: 



1 f"^ xe^dx 

1 xe^lg^ - e^^lg^ _ 1 e°^/2(aA - 1) + e^''^/^ 
A ~ 2A sinh(aA/2) 

In most general cases A « 10^^ is a big number for ceramic materials, so the value of Qx can be calculated as A tends 
to infinity. 

1 e''^/2(aA- 1) + e-"^/2 ^ ^ 

lim 77": . , , , = a (18) 

x->oo 2A sinh(aA/2) ^ ^ 

Thus we see that in limes for small misfit angles and the condition that cos tp 3> sin ip we have uniform thickness of 
the intergranular films. This is in good agreement with the experiments. Lets calculate the special case where A = 0, 
this the case when angle of misfit is 6* = 0. 

1 e''^/2(aA- 1) + e-'»^/2 

lim 77T = a/2 19 

A->o 2A sinh aA/2) ' ^ ' 
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Previous calculation can be used together with Eq.ljl]) to calculate the thickness (d) — of intergranular film in 

polycrystalline material. From this consideration we conclude that average thickness over wide range of crystallo- 
graphic orientations is almost constant and it reduces to half this value only in the case where the misfit angle is 
zero. 



V. KINETICS OF THE INTERGRANULAR FILM 



The kinetics of interface in one dimensional case is allready given by Mori et al.— Solution of Eq. ([2]) is given as: 

{p/a)v + 1 

Where v — Vr/a is non-dimensional velocity of the interface and V is real velocity of the interface. Velocity diverges 
at point where /3 — *■ a this is non-physical solution but the validity of Time-Dependent-Ginsburg-Landau equation 
(TDGL) is questionable in the range of very high velocity of the interface. We can give several limiting cases for 
Eq. (1201): 

V > a/T {(3 / a): 

..-l/21ni^M^ 
' 2 

-V > a/T {(3 \ -a) : 

1/91 1 + ("//^) 
V ~ 1/2 In 

We conclude with application of evolution equation to phase separation in polymer mixtures. This is usually 
described with Flory-Huggins free energy given as: 

F 

— = (t)A\n(t)A + (j)Bln(j)B + Nx(l)A(t)B (21) 

Where N is degree of polymerization, x interaction parameter per monomer (Flory-Huggins parameter), (j)A and (j)B 
is probability of site in space is occupied with molecule of A or molecule of B species, respectively. Lets take a 
characteristic length of polymer chain to be is given hy w = a\/N where a is size of monomer and N characteristic 
loop size of random walk of polymer A protruding in space of polymer B. Then total energy coming from interaction 
between A and B segments of polymer \s L = xNkT so a = = xN. If we put this result inside Eq. ^ we have: 



^^'2w _ 2a^Nioo,p _ 2a ^^^^ 
a xNioop Xy/Nioop 

For equilibrium the interaction energy of monomer will be of order kT thus: X-^ioop ~ 1- Giving finally in equilibrium 
case: 

Which is of course not a surprise as Temkin's model is equivalent to Flory-Huggins theory of polymer mixtures. 
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VI. APPENDIX 



SI _ dl _d_ I df_\ 
^T- = {{ax(l>x + ay(j)y + az(j)z) ln(aa:0x + ay<Py + az(l>z)) (25) 



Then 



d_ / 5/ 
dx \d(f>x 

d ( 2^1 : 1 . 1 



ln(aa;(?!)a; + a^c/ij^ + az^z) + aa:'/'^;— =77 + ay(t)y——ax + a 



OX \ av(p av(j) avcp 



a. 



+ ay(f)yx + az4>zx) + al^xx-^ - alcf) + ay (pyx + az4>zx) 



) + Uzf^zx ^i ^x 

a\(p [avcp)'' a\(p 

- az4>z ^^Y^^2 iP'x4>xx + ay (pyx + az4>zx) 

{P'x4*xx + ay(j)yx + az4>zx) - , ^ <\>{a-x4>xx + ay4>yx + a-z^zx) 



v-fi\'^x-rxx I ^yyyx ' ^z^zxf / t-tjL\0 

aV0 {a\/(py 

a, 

{flx^Pxx 



+ -z^iaxfpxx + aycpyx + az<t>zx) 



aV(f) 

Similarly, for components y and z we have : 



-( 




ay 


dy \ 


^d(j)y) 


aV(j) 


d 


(df\ 


1 az 


dz 







Then the sum of these terms is: 



dx \d(j)xj dy \d(j)yj dz xdcpz 

^^(«x^^xx + al(pyy + al4'zz + ^ttxttyCpxy + ^axaz<Pxz + ^ayaz4>yz) 

aV(t> 



(26) 



(27) 
(28) 

(29) 



A. One-dimensional solution to evolution equation 



We start from partial differential equation: 



with change of variables: u{(j)) = a-^ we have: 



du(<j)) d ( d'^cpdx d'^cj) A 



d(t> \ dxj dx^ d(j) dx^ dx 



|jj = — and we transform the equation in : 

du 



vu — a(f) + p = 



uU) = Ce-""^ - -(- - ^) - - 

V V V 



where 



Final solution is: tanht; = (^i3/^)v+i '^liich for /3 = gives tanht; = —v. 



B. One-dimensional solution of evolution equation with harmonic oscillations 

^=;3-a</, + a(l + cos(a;i))0/|! 
Lets solve it by separation of variables: 

(j){x,t) = X{x)T{t) 

Where X(x) is only function of x and T(t) is only function of t. Then, derivatives are: 



dx"^ dx^ 
d(i) ^, , dX 

d<i> dT 



Substituting these derivatives in differential equation: 



Dividing this differential equation with X we get: 

1 , , ,d^X,dX , dT 

-a(l + co,sM)^/--.T + /3A=- 

If search for solution of where there is no driving force for transition so we have /3 = 0: 

1 d'^X ,dX , ^ dT^ 1 

-/^ - {<^T+— )—-——- = 



X dx dx dt a(l + cos(a;t) 
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This equation is separable, and we get two ordinary differential equations: 

(fx ^ ^dX 



dx'^ dx 

Where m is integer: 



- 2mX— = (42) 



dT 

— +aT- 2ma{cos{(jjt) + 1) = (43) 
The solution of second differential equation is simply: 

T{t) = C exp(-at) + — — ——-a{a'^ + J^) (44) 

The solution to first differential equation can be found by putting solution X(x) back in the equation: 

X{x) = - tanh(m(a; - 6)) (45) 

So final solution is: 

2am(a^ + + cos{ut) + auj sva.{ojt)) 



(p{x, t) — — tanh(TO(a; — 6)) 



C. Calculation of Read-Schokley integral 



(46) 



Calculation of the integral: 



= ^Jdw- = ^''^ 

J cos(^exp [Ga(cos0 + sin(/))/(47r(l - (T))d{A - ln6)/kT] sin 



/ cxp [Ga(cos0 + sin (/))/(47r(l -tT))6'(^- In 6l)/A:T] sin 

We will evaluate this integral with method given in Massida. All constant factors can be extracted from previous 
integral and then we need to solve general integral of type: 

/■27r 

/ e''"^exp{cos(f) + sin (l))d<l) (48) 
Jo 

Using relation between exponential and Bessel function: 

^ k—oo 

exp(|(t + -)) = (49) 

k— — oo 

gcos = ^ (1) = J2 e'^'^h (1) (50) 

fe k 

gSin0 ^ ^(-i)''e^'"^4(l) (51) 



Thus integral can be written as: 



/ e'""^ V(-i)''e^'''^4(l) Ve^*='^4(l)d6l (52) 
Jo h ^ 



Using = Im{x) we can transform integral to 2n J2k Ik{^)Ik+i (~*)'^ using addition theorem for Bessel functions: 

^(_l)ne-*J„(^)/„+,(Z) = I,{Lo)e'-^ (53) 
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Where w = (z^ + Z'^ — 2zZcos(\)Y^'^ and is defined by Z — zcosc/) = u;cos-0 and zsinc/) = wsin?/". The addition 
theorem can then be applied to k summation in sum and we get solution as: 



27r^4(l)/fe+i(l)(-i)'= = /i(v/2(l - cos,/,))e^™(V2(i— «/2) (54) 



fe 

Final solution to integral is: 



cos(x) exp [cos(a;) + sin(x)] dx — 2-Kl\{\f2)—^ (55) 







Similarly for: 



1/2/ sin(2x)exp[cos(a;)+sin(a;)](ia; = /„(w)e'^"^ = /2(V2)27r (56) 







D. Two-dimensional solution of evolution equation 

We start from partial differential equation: 



dr 

Using transformation from Nerney et ali^ we can write previous equation: 

^ - a0 + /3 = -w6'(0) (58) 
d(p 

Where v — non-dimensional velocity, solution to eqution is: 

e = Ce"^ + -0-/3/f -a/w^ (59) 



V 



Where C=\/(-% + -) — ( — For equilibrum conditions is equation reduced to: 



Solution of partial differential equation is: 

r dr a 



^ ^ (r0) = -e*« (61) 



Determination of constant of intergration C: 



j 1 when r = 0, 

(r) = < . , (62) 
— 1 wlien r = 00 



Giving C — Q and finally the solution to differential equation: 

• In 



0W = (63) 



